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ABSTRACT 

We use a combination of a cosmological A-body simulation of the concordance A cold 
dark matter (ACDM) paradigm and a semi-analytic model of galaxy formation to 
investigate the spin development of central supermassive black holes (BHs) and its 
relation to the BH host galaxy properties. In order to compute BH spins, we use the 
a model of Shakura & Sunyaev and consider the King et al. warped disc alignment 
criterion. The orientation of the accretion disc is inferred from the angular momentum 
of the source of accreted material, which bears a close relationship to the large-scale 
structure in the simulation. We find that the final BH spin depends almost exclusively 
on the accretion history and only weakly on the warped disc alignment. The main 
mechanisms of BH spin-up are found to be gas cooling processes and disc instabil- 
ities, a result that is only partially compatible with Monte-Carlo models where the 
main spin-up mechanisms are major mergers and disc instabilities; the latter results 
are reproduced when implementing randomly oriented accretion discs in our model. 
Regarding the BH population, we find that more massive BHs, which are hosted by 
massive ellipticals, have higher spin values than less-massive BHs, hosted by spiral 
galaxies. We analyse whether gas accretion rates and BH spins can be used as trac- 
ers of the radio loudness of active galactic nuclei (AGN). We find that the current 
observational indications of an increasing trend of radio-loud AGN fractions with stel- 
lar and BH mass can be easily obtained when placing lower limits on the BH spin, 
with a minimum influence from limits on the accretion rates; a model with random 
accretion disc orientations is unable to reproduce this trend. Our results favour a sce- 
nario where the BH spin is a key parameter to separate the radio-loud and radio-quiet 
galaxy populations. 
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1 INTRODUCTION 

Galaxies with Active Galactic Nuclei (AGN) have become 
the centre of attention in extragalactic studies due to their 
star role in the galaxy formation scenario favoured today. 
They are a key ingredient in explaining several observational 
statistics on galaxy star formation rat es, lum inosities and 
colours (see for instance, iBo wcr ct al. 1 20061 . ICroton et al 



20061 . ICattaneo et al . 2006, Si iacki et al1l2007l. iMarulli et al 



20081 . lLagos. Cora fc Padillal 120081 . ISomerville et al.l [2008). 
However, a number of phenomena related to AGN are still 
unclear, including the origin of the AGN dichotomy into "ra- 
dio loud" (RL) and "radio quiet" (RQ) types. RL AGN are 
characterised by the presence of jets of relativistic plasma 
and/or associated large-scale radio lobes. Relativistic jets 
emitting synchrotron radiation in the radio band are gener- 



ated by the collimation of outflows from the innermost re- 
gions of accretion discs around black holes (BHs), carrying 
away a large fraction of the available accretion power. How- 
ever, the jet production mech anism is not fully under stood. 
The double-lobed morphology Fa naroff fc Ril ey (1974) Class 
II (FR II) sources are the most common type observed in 
broad line emitting RL quasars. On the other hand, the 
lower luminosity Fanaroff Class I objects (FRI) have very 
weak radio emitting ejecta and weak or no emission lines. 
The correlation between radio luminosity and line luminos- 
ity in [OIII] A5007 is sim ilar for both RL and RQ AGN 
l|Xu. Livio fc Baumll 19991 ). 

Different criteria for the definition of radio loudness 
have been used, the most com mon being the ratio be tween 
radio-to-optical l uminosity. Ti |Kellermann et al.lll989l ). and 
the radio power (|Miller et al.lll990h . A third criterion based 
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on th e classical radio mo rphology has also been consid- 
ered dSulentic et al.l l2003h. From the combination of the 
two latter criteria. Izamfir. Sulentic fc Marzianil (|200ST l find 
that the boundary between RL and RQ AGN is given by 
log Li. 4 ghz = 31.6 ergs - Hz -1 , which corresponds to the 
luminosity of the weakest FRII sources. Using this luminos- 
ity threshold, they obtain a fraction of ~ 5 — 9 per cent of RL 
AGN relative to the total population of AGN. The presence 
of a physical bimodality in the properties of AGN is still a 
matter of debate, and it should be borne in mind that ob- 
serv ational selections might actually induce t his dichotomy 
(e.g. ICirasuolo et~ai1l2003l . IWhite et al]l2007l ). 

The host galaxies of both RL and RQ QSOs with 
nuclear luminosities My < —24 a re massive ellipticals 
ijDunlop et al.ll200llFlovd et al.ll2004h ; disc-dominated host 
galaxies become increasingly rare as the nuclear power in- 
creases. Sikora, Stawarz & Lasota (2007) demonstrate that 
radio-selected AGN hosted by giant ellipticals can be up to 
a thousand times brighter in radio frequencies than AGN 
hosted by disc galaxies, thus generating two sequences in 
the radio versus optical luminosity plane. These RL and RQ 
sequences remain present when considering the dependence 
of the radio loudness 1Z with the Eddington ratio, where 
the latter is defined as the bolometric luminosity (estimated 
from optical data) in units of the Eddington luminosity. The 
two sequences show an increase of radio loudness with de- 
creasing Eddington ratio. This bimodality breaks down at 
hi gh accretio n rates. This behaviour was originally noticed 
bvlHo1 l|2002h . although only a single sequence was visible as 
a consequence of the reduced sample of AGN considered. 

The link between jet production and gas accretion rate 
onto the BH has been widely investigated. The inner part of 
the accretion flow is believed to be responsible for the ob- 
served hard X-ray emission produced by Compton scattering 
or Bremsstrahlung radiation. Thus, observations of the com- 
pact emission in the X-ray and radio bands help understand 
the disc-jet connection i n both stellar and supermassive BHs 
dUlvestad fc Hoi l200ll. iMerloni. Heinz fc di Matted 12003! 



Falke. Kording fc Markofj |2004 iKording. Jester fc Fender! 
20061 ). under the assumption that the jet formation process 



is qualitatively similar on both scales. Observations of black 
hole binary systems during the low/hard X-ray state reveal 
a tight correlation between th e radio and X-ray luminosity 
jGallo. Fender fc Poolevl 1 2003 ) , supporting the dependence 
of the radio power on the accretion rate. This "acretion 
paradigm" postulates that radio loudness is entirely related 
to the states of accretion discs. However, it cannot explain 
the two parallel sequences AGN occupy in the plane de- 
fined by rad i o lou dness and Eddington ratio, as found by 
ISikora et all (|2007l ). Hence, it is possible that additional BH 
properties could be playing an important role in driving the 

production of powerful relativistic radio jets. 

The RL and RQ sequences found bv lSikora et~ail (|2007l ) 
are characterised by AGN containing BHs with masses 
> 10 s M Q and < 3 x 1O 7 M , respectively. Nonetheless, 
the BH mass is not capable of separating completely these 
two sequences, since objects powered by equally massive 
BHs ca n di f fer in radio luminos i ty by ~ 4 orders of mag - 
nitude jHcl l2002l . iDunlop et ail I2003L ISikora et al] 120071 ). 
iBlandfordl l)l990h introduces a new parameter suggesting 
that the efficiency of jet production is determined by the 
dimensionless BH spin, according to the assumption that 



relativistic jets are powered by rotating BHs through the 
Blandford-Znajek mechanism. In an attempt t o explain this 
radio power dichotomy at fixed accretion rates, ISikora et al.l 
(2007) propose a revised version of the "spin paradigm" 
which suggests that giant elliptical galaxies host, on av- 
erage, BHs with larger spins than those hosted by spiral 
galaxies. These new insights inferred from observational ev- 
idences motivated the study of the influence from BH merg- 
ers and small accretion episodes on the evolution of BH 
spins and their relation to the morpho l ogy of the host g alaxy 



and their relation to tne morphology 01 the nost galaxy 
Volonteri. Sikora fc Lasota! 120071 ; iNemmen et alj|2007l ; 



iBerti fc Volonteril l200Sf ) . These theoretical works indicate 
that it is feasible that the observed morphology-related bi- 
modality of AGN radio loudness arises from the dichotomy 
shown by the distribution of BH spins. 

In this paper, we investigate the latter possibility by 
using the SAG (acronym for 'Semi- Analytic Galaxies') semi- 
analytic model of galaxy formation by Lagos, Cora & Padilla 
(2008, hereafter LCP08), combined with the outputs of a 
cosmological simulation based on the ACDM cosmology. The 
advantage of this approach relies in that we are able to use 
realistic BH growth histories directly linked to the evolution 
of the star formation rate of their host galaxies. The SAG 
model also allows us to evaluate the influence of the different 
mechanisms that contribute to the accretion of material onto 
the BHs, which drive both the BH mass growth and the BH 
spin development. 

This work is organised as follows. In Section [2] we 
briefly summarise the mechanisms that contribute to the 
BH growth included in SAG, and describe in detail the im- 
plementation of the BH spin model. In Section [3] we discuss 
the connection between the BH spin distributions and the 
characteristics of the BH host galaxies. Section [3] shows the 
results on the final BH spin and the way in which it is ac- 
quired. In Section [5j we analyse the nature of AGN radio 
loudness and compare our results with available observa- 
tional data. Finally, the main conclusions obtained in this 
work are summarised in Section [6] 



2 COMPUTING THE BLACK HOLE SPIN 

We study the relationship between the BH spin and the 
properties of the BH host galaxy using a combination of a 
cosmological TV-Body simulation of the concordance ACDM 
universe and the SAG semi-analytic model of galaxy for- 
mation described in LCP08. This model considers radiative 
cooling of hot gas, star formation, galaxy mergers, disc insta- 
bilities, feedback from supernova explosions, BH growth and 
AGN feedback produced during accretion onto BHs driven 
by gas cooling processes. 

In the following subsections, we briefly describe the TV- 
Body cosmological simulation and the BH growth model 
used in SAG, relevant for the present study. We then ex- 
plain in detail the implementation of the model to compute 
the BH spin. 



2.1 ACDM Cosmological simulation 

The cosmological simulation is the same as the one used by 
LCP08, based on the concordance ACDM cosmology. It con- 
siders a periodic box of 60 Mpc containing 16, 777, 216 
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dark matter particles with mass 1.001 x 1O 9 /i _1 M . It has 
more than 54, 000 dark matter haloes with masses of up to 
5.36 x 10 14 M©. The simulation p arameters are consis - 
tent with the results of WMAP data (|Spergel et al.ll2003h . 
with f2 m = fiDM + fibaryons = 0-28 (with a baryon fraction of 
0.16), Q,a ~ 0.72 and as = 0.9. The Hubble constant is set 
to Ho — 100 /iMpc -1 , with h — 0.72, and the gravitational 
softening length is e = 3.0/i -1 Kpc. The simulation starts 
from a redshift z = 48 and has been run usin g the public 
version of the GADGET-2 code (|Sprmgej|20"05T ). 

2.2 Black hole growth 

There are three distinct modes of BH growth in SAG. The 
first one, referred to as the "QSO mode", is associated to 
starbursts in galaxies. These events can be driven by galaxy 
minor or major mergers and the collapse of unstable discs. 
The second one, referred to as the "radio mode", is associ- 
ated to gas cooling processes that produce star formation in 
a more quiescent way. The third one is BH mergers which 
occur shortly after the merger of the parent galaxies. In our 
model, AGN feedback occurs only during the latter mode 
which produces low accretion rates, in agreement with ob- 
serva tions (e.g. lHoll2002l ; iDonahue et~aill2005l ; ISikora et all 
120071 ). The detailed assumptions for each mechanism are 
given in LCP08. Here, we briefly describe the main formulae 
involved. 

During galaxy mergers, the BH growth depends on the 
ratio bet ween total galaxy ma ss and the amount of available 
cold gas jCroton et alJ :2006, LCP08) as 



Mb 



/bh M" 
At M cen 



McoldGa 



1 + (200kms-VKir 



(1) 



where McoidGas is the cold gas mass of the system, consist- 
ing of the central and satellite galaxies, with masses (in- 
cluding stellar and cold gas mass) M sat and M ccntral , re- 
spectively; At is the timescale between consecutive SAG 
timesteps. Black holes are assumed to grow with an effi- 
ciency /bh = 0.015, set to fit the Mbh-MbuI^c relation. 
For disc instabilities, the ratio M sat /M contral is replaced by 
unity, since this process depends on the properties of a single 
galaxy. It is important to note that there are no BH seeds in 
our model and, therefore, the birth of central supermassive 
BHs is triggered when galaxies undergo their first starburst 
at some point in their evolution. 

The accretion rate during the "radio mode" is given by 

3 

(2) 
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,V MbH /hot 

MbH = KAGN 1 „„ ; X — — X , 

10 8 M Q 0.1 \200kms -1 

where Mbh is the BH mass, and /hot is the fraction of the 
total halo mass in the form of hot gas, /hot = wiHotGas/M v i r , 
Mvir being the virial mass of the host halo which has a 
virial velocity Kd r ; kagn is a free parameter set to fit the 
galaxy and QSO luminosity functions with a value 2.5 x 
lO^Moyr" 1 . 

Irrespective of the accretion mode, the BH produces 
a luminosity that corresponds to the total release of en- 
ergy from the disc formed by gas accretion. This lumi- 
nosity is given by Lbh = V Mbhc 2 , where c is the speed 
of light, and r\ — 0.1 is the standard efficiency of energy 
production that occurs in the vicinity of the event horizon 
jShakura fc Sunvaevlll973l SS73 hereafter). 



2.3 Black hole spin 

Our implementation of a BH spin model in SAG follows 
the same methodology adopted for all the galaxy proper- 
ties considered in the semi-analytic model, which consists in 
solving a set of linearised differential equations. The BH spin 
in each galaxy is followed through the evolutionary history 
of the simulated universe, allowing us to analyse the way 
in which the BH spin is acquired for every host galaxy and 
to study directly the role and influence of the hierarchical 
build-up of structure in the development of this BH prop- 
erty. During galaxy mergers, the two central supermassive 
BHs are also merged; we simply consider that the mass of 
the final BH is the sum of the merging BH masses; the re- 
sulting spin is calcu lated using the semi-a nalytic fitting for- 
mulae presented bv lRezzolla et ail (|2008l ) to fit a complete 
compilation of numerical results on BH mergers. The other 
processes involved in the BH growth feed baryons through 
the formation of accretion discs. The treatment of the an- 
gular momentum transfer from the accretion disc to the BH 
requires the following analysis. 

The dimensionless BH spin is defined as a = 
./bh/Jmax = cJbh/GA/| h , where Jbh is the angular mo- 
mentum of the BH. In general, J refers to the absolute 
value of the angular momentum, \J\. We consider the pres- 
ence of a warped disc characterised by the angular momen- 
tum of the accretion disc at the warp radius, Jd = Jd(Rw) 
l|Volonteri. Sikora fc Lasotall2007l '). In this case, the amount 
J d/ 'Jbh can be used to define the alignment betwee n the 
spinning BH and the accretion disc (|King et al.ll2005l . here- 
after K05). Angular momentum vectors are co-aligned when 
cos((f>) > — 0.5Jd/ Jbu, being <j> the angle between the disc 
angular momentum and the BH spin; otherwise, the accre- 
tion disc is counter-aligned with the BH spin. We use this 
criterion to evaluate the BH spin-up or down in the presence 
of a warped disc. 

SAG does not assume an initial mass for BHs, or 
BH seeds; instead, the first starburst in a galaxy, that 
occurs during galaxy mergers or disc instabilities, cre- 
ates the central BH. These processes would therefore be 
in line with an origin of super-massive BHs as rem- 
nants of the first events of stel lar explosions early in 
the hi story of the Universe (e.g. Bernadetta fc Volonteril 



|200S|. lOmukai. Schneider fc Haimanl l2008llShapirol 120051 
lElmegreen. Bournaud fc Elmegreenl 120081 ). Given this pos- 
sible stellar origin of the BHs, it is reasonable to expect 
their initial spin to be non-zero. However, for the sake of 
completeness we also study the effect of initially spinless 
BHs, which could b e formed by rapid collapse of inf alling 
gas (see for instance iBegelman. Volonteri fc Reesll2006h . In 
this latter case, BHs acquire spin for the first time during 
their first accretion episode or merger with other black hole. 
In the case of an accretion even t, the res u lting spin is given 
by the formalism developed by iBardeenl (|l970h . When two 
spinless BHs mer g e, we use the analytic formula given by 
iBerti fc Volonteril £2008); when at least one merging BH has 
a non-zero spin, the resulting spin is calcul ated using the 
semi-a nalytic fitting formulae presented by iRezzolla et al.l 
(2008) to fit a complete compilation of numerical results on 
BH mergers. 

In the case of BHs with non-zero spin, we calculate Jd 
taking into account the mass of the accretion disc, Md, and 
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the remaining properties related to the BH, A/bh, Mbh and 
Lbh; all these quantities are provided by the SAG model. 
The ratio between the angular momentum of the accretion 
disc and the BH spin is expressed as 
1/2 

(3) 



Jd 
Jbh 



M d 



aM BH 



where = MbhA^ The accretion rate Mbh comes either 
from equation [T]or [2] Accretion discs are completely con- 
sumed in the time interval At. The ratio R w /R s is the warp 
radius expressed in terms of the Schwarzschild radius, R s ; 
in a SS73 middle-region disc, this ratio can be written as 
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(e.g. K05. IVolonteri et al.l 120051 ). where /ecw = iBH/^Edd, 
being LEdd the Eddington luminosity, vi the warp propaga- 
tion viscosity, and v\ the accretion driving viscosity (1/2 can 
be different from v\). The parameter a, introduced by SS73, 
takes into account the efficiency of the mechanism of angular 
momentum transport. A thin accretion disc is characterised 
by H/R < a < 1, where H and R are the disc thickness 
and radius, r espectively. In this case, it c an be shown that 
U2/U1 ~ a 2 (|Papaloizou fc Pringle|[i"983h . Nonetheless, at 
high accretion rates (i.e. hot /thi ck accretion discs) there are 
indic ations that a > 0.1 (e.g. iDubus. Hameurv fc Lasotal 
120011 ). Based on these theoretical and observational indica- 
tions, we adopt simple criteria to distinguish between thin 
and thick accretion discs and, therefore, to choose an appro- 
priate value for a: 

• Thick accretion disc: if the accretion occurs at super- 
Eddington rates, that is Lsn/LEdd > 1, the accretion pro- 
ceeds vi a a radiation pressure supported thick disc with H ~ 
R (e .g. IVolonteri etail 120051 ; iBegel man. Volonteri fc Reesl 



I2006T ) , which imp lies v\ V2 (e.g. Kumar fc Pringlel 19851 ) 



and a > 0.1 (e.g. lDubus. Hameurv fc Lasotall200ll ). 

• Thin accretion disc: if the accretion occurs in the regime 
of Z/BH/iEdd < 1, it is assumed that H/R ~ a C 1. 

In order to test the reliability of our results, we adopt 
different values of the parameter a, ranging from ~ 0.1 to 
~ 0.5, for the thick disc, and from ~ 10~ 4 to ~ 0.05, for 
the thin disc. It is important to remark that the distinc- 
tion between thick and thin discs has a direct consequence 
on the amount of angular momentum per unit mass trans- 
ferred from the accretion disc to the BH, for the same disc 
mass. This arises from the dependence of equation [4] on the 
parameter a, and the viscosities v\ and V2- This is a simpli- 
fied approach to more sophisticated model s which emphasize 
the importance of the disc thickness (e.g. lHawlev fc Krolikl 
l200fj ). 

There are a number of works that study in de- 
tail the possibility of alignments between galactic- 
scale structures and the ce ntral BH-disc system. 
iBogdanovic. Reynolds fc Millerl (|2007l ) find indications 
of co-alignment between the BH spin a nd large-scale gas 
flows coming from galaxy mergers. Also, iLindner fc Fragile] 
(|2007i ) show that the jet direction is set by the BH 
angular momentum. These results indicate a possible 
co-alignment between the jet direction and the angular 
momentum of the accretion disc of the host galaxy, also 
supported by observations of high redshift radio galaxies 



^Chambers. Milev fc van Breugell Il987l; lLacv et al.l ll99S 



llnskip et all 120051 ). Moreover. iBerti fc Volonteril (|200S ) 
show that high spin values are efficiently achieved by a 
model with non-chaotic accretion episodes, that could 
explain t he very large spin a ~ 0.99 of the case of MCG- 
06-30-15 ((1 rcnncman fe Reynolds! 120061 ). These arguments 
make it relevant to consider the possibility that the direction 
of the angular momentum of the accretion disc is related to 
the physical process that triggers the gas accretion. 

We test this possibility using the information provided 
by our model as follows. When the accretion is driven by disc 
instabilities, we assume that the direction of the accretion 
disc can be inferred from the angular momentum LGai,disc 
of the galaxy disc, 



Jd = Jd ■ ^Gal.dis 



(•») 



On the other hand, if the accretion is driven by gas cooling, 
we assume that the direction will be that of the angular 
momentum of the host halo, 



Jd = Jd ■ Lf 



(6) 



If the gas accretion is driven by galaxy mergers, the direction 
is taken from the gas mass- weighted disc angular momentum 
contributed by each galaxy, 



Jd — Jd 
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(7) 



where Mg^Q^ s and M^ MG ^ S are the cold gas mass of the 
central and satellite galaxies, respectively. Having the in- 
formation on Jd, we determine the co- (counter-) alignment 
according to the K05 criteria to finally sum (subtract) Jd to 
Jbh. 

Since a ACDM universe is characterised by strong co- 
alignments between the angular momentum of dark-matter 
haloes an d their surrounding structures on scales of sev- 
eral Mpc l|Paz. Stasvszvn fc Padillal l2008h . there could be 
some level of alignment between the BH spin and the large- 
scale structure. In order to explore this, the initial angular 
momentum of the gaseous disc of our model galaxies are 
inferred from that of their host dark-matter haloes. Then, 
the angular momentum of model galaxies is followed and 
updated at every star formation episode, which can intro- 
duce some degree of misalignment with that of the hosting 
haloes. In the case of gas cooling processes, we consider that 
the cooled gas added to the g alaxy disc drags the angular 
momentum of the host halo (|Cole et alj |2000| ) , modifying 
the pre-existing angular momentum of the galaxy disc. Dur- 
ing mergers and disc instabilities, the contributions from 
the participating gas and stellar components are taken into 
account to obtain the resulting galaxy angular momenta. 
We also consider ran d omly oriented accretion events (e . g . , 
iKing fc Pringlei 120071 ; IVolonteri. Sikora fc Lasotal 120071 ) in 
order to study the main consequences arising from these 
two possible scenarios. 

We also analyse the results obtained from a model 
where the timescale for disc warping is much longer than 
for gas accretion. In this case, we ignore the K05 alignment 
criterion and the angular momentum of the disc is simply 
vectorially added to the BH spin. In all cases, the final BH 
spin direction will be that of the total angular momentum of 
the system (BH and accretion disc) which is assumed to be 
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conserved. This assumption is not fully accurate since the 
torque exerted by the accretion disc acting over the total 
angular momentum of the system may be important (e.g. 
IVolonteri et~alll2005h . 

Since we are not considering angular momentum losses 
by torques, our estimates of the values of BH spins are upper 
limits. In case that the spin takes unphysical values (a > 1), 
we simply set them to a = 1. However, when the resulting 
spin values are more than half a decade larger than unity, 
we also reduce the gas accretion onto the BH consistently 
with this new upper limit. We find that adopting values of 
a — 0.5 for thick discs and a = 0.05 for thin discs, only 15 
per cent of all accretion events in our model produces values 
of a > 1 (only « 0.1 per cent yields a > 3). Lower values of a 
produce higher percentages of unphysical spin values, as can 
be inferred from equations [3] and [4] Therefore, our choice of 
a produces a minimum bias on the resulting spin distribu- 
tions (a high concentration at a~ 1 would spuriously hide 
information on the high tail of the BH spin distribution). 

One of the main uncertainties in our spin model lies 
on the initial value of the BH spin when considering their 
possible stellar origin, for which we test two possibilities, (i) 
a constant initial value (from Sinitial = 10 -4 to O.fQ), and 
(ii) an initial value proportional to the BH mass. As we will 
show in Section[3l the choice of the initial spin value does not 
affect our results significantly. Moreover, the general conclu- 
sions of this work are preserved even for the case of spinless 
initial BHs. The most sensitive population to the inclusion 
of an initial BH spin are low-mass BHs (i.e. Mbh < 1O 6 M0) 
since most of them have very few accretion episodes and 
practically no merger events during their entire lifetimes; as 
a note of caution, we remind the reader that a null initial 
spin is not likely to occur in a BH with a stellar origin. 

In the case when a BH is assigned an initial spin (ap- 
part from a merger event) the spin direction will depend on 
the process that triggers the BH birth. Thus, the direction 
is given by (i) the angular momentum of the galaxy disc 
(similar to Eq. [BJ, in the case of disc instabilities, and (ii) 
the angular momentum of the galaxy disc weighted by the 
gas mass of each merging galaxy (similar to Eq. [7}, in the 
case of galaxy mergers. The final value of the spin will be 
only a function of the accretion and assembly history of the 
BH. 

This model provides a complete characterisation of a 
supermassive BH, giving information on both its mass and 
spin. These quantities strongly depend on the mass of cold 
gas accreted onto the BH. The great advantage of our model 
is that accretion discs are consistently given by our semian- 
alytic model within a universe where structures grow hier- 
archically, instead of being obtained from the assumption 
of specific distrib utions of Eddington accretion ratios as in 
other works (e.g. ISikora et al.ll2007h . Indeed, the redshift 
evolution of the accretion rate density onto a BH found 
by LCP08 (see their Fig. 2) is consistent w i th th e empir- 
ical estimations given by iMerloni fc Heinj l|200ct ). More- 
over, the evolution of the QSO luminosity function obtained 
from SAG is in good agreement wit h observational results 
jWolf et alJl2003 ; ICroom et alj|2004 ). Consequently, the re- 



Thesc initial spin value s are consisten t with results from hydro- 
dynamic simulations (see [Shapiro 200^ and references therein). 




logio(M Bu ig 8 [MJ) 

Figure 1. Distributions of z = bulge masses for different ranges 
of final BH spin values, a< 0.5 (dashed lines) and a^ 0.5 (solid 
lines). Black lines correspond to a model where accretion discs 
orientations are inferred from galaxy discs according to the pro- 
cess that triggers the accretion event; grey lines denote random 
accretion discs orientations. Both models take into account the 
K05 alignment criterion and use a constant initial black hole spin 
value a; n j t ; a i = 0.01. 



suiting distributions of accretion disc masses and Edding- 
ton rates are consistent with what was adopted in previous 
Monte-Carlo models. 

2.4 Accretion disc orientation: inferred or 
random? 

An important improvement presented by this model resides 
in the use of the known, vectorial angular momenta of the 
components of model galaxies to infer the direction of the 
angular momentum of the accretion disc. This depends on 
the physical process that triggers the accretion event, that 
is, gas cooling, disc instabilities, or galaxy mergers, as de- 
tailed above. Here we describe the main advantages of this 
approach. 

Fig. [T]shows the distributions of bulge masses for galax- 
ies hosting central BHs with low and high spin values at 
z = (dashed and solid lines, respectively), for a model 
where the angular momentum of the accretion disc depends 
on the process responsible for the gas accretion onto the 
central BH (Eq. [6] and (model A, black lines), and a 
model that uses randomly oriented accretion discs (model 
B, grey lines). Both models take into account the K05 align- 
ment criterion and use a constant initial black hole spin 
^initial = 0.01. As can be seen, in model A, the fraction of 
galaxies with bulges of masses Meuigc >3x 1O 9 M0 hosting 
low spin BHs decreases faster than in model B, while both 
models show similar fractions of galaxies for a given bulge 
mass that host BHs with high spin values. Results from both 
models indicate that massive bulges host exclusively high 
spin BHs. Taking i nto account tha t the most massive galax- 
ies are ellipticals (jConselicel 12006). we find that model A, 
which shows a slightly clearer connection between massive 
galaxies and high BH spin values, supports the conjecture 
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Figure 2. History of average BH growth for elliptical and spi- 
ral galaxies (top and bottom panels, respectively). Different lines 
show different average final BH masses (each mass bin contains 
equal number of BHs). Elliptical galaxies are selected according 
to MBuigo/MgtciiarTotal > 0.95, while spirals satisfy the condition 

< M B ul g e/MstellarTotal < 0.95. 



proposed bv lSikora et al.l (|2007l ) that the radio loudness bi- 
modality can be explained by the morphology-related bi- 
modality of the BH spin distribution. Sections [3] and [3] are 
devoted to explore the predictions of this model. 



3 CONNECTION BETWEEN BLACK HOLE 
SPIN DISTRIBUTIONS AND HOST 
GALAXY PROPERTIES 

The main purpose of this work is to find a link between the 
properties of BHs and those of their host galaxies. In order 
to test the reliability of the possible relations obtained, we 
evaluate their dependence on different physical assumptions 
included in our BH spin model, such as the initial value of 
the BH spin and the alignment of the accretion disc. 

Fig. shows average growth tracks for BHs hosted by 
elliptical (top) and spiral galaxies (bottom). In both cases, 
we show BHs with masses Mbh > 10 Mq , separated in 
three mass bins, each containing the same number of BHs. 
The morphological criterion to distinguish between ellipti- 
cal and spirals considers the fraction MBuige/MsteiiarTotai 
between the bulge mas and the total stellar mass of the 
galaxy. We consider as ellipticals galaxies, those with r = 
Msuige/A/stciiarTotai > r t hresh = 0.95, chosen in order to 
recover observed morphological distributions (see LCP08). 
These tracks indicate that, at all times, ellipticals host more 
massive BHs than spirals. Besides, the latter tend to stop 
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Figure 3. Distributions of black hole spins at z = hosted by el- 
liptical (left panels) and spiral galaxies (right panels) from model 
A, in which the direction of the accretion disc angular momen- 
tum is inferred from the process responsible for the gas accre- 
tion onto the central BH. The distributions are obtained from a 
constant initial BH spin a; n ; t ; a i = 0.01 and the K05 alignment 
criterion (upper panels), from a initial = M B h mitial/10 7 MQ and 
K05 (middle panels), and from a; n ; t iai = 0.01 and by adding the 
angular momentum of the accretion disc to the BH spin without 
considering the K05 criterion (lower panels). The vertical solid 
lines represent the average values of BH spin obtained from the 
histograms. For reference, the average values of BH spin obtained 
from randomly oriented accretion discs (model B) are shown as 
vertical dashed lines. 



the growth of their central BH at earlier times. These results 
do not depend on the details of the BH spin model. The dif- 
ferences in final mass and assembly history of the BHs can 
produce different final BH spins, a possibility that we now 
turn to. 

Fig. |3] shows the BH spin distribution for BHs hosted 
by elliptical (left panels) and spiral (right panels) galaxies 
for three different versions of model A, (i) using a con- 
stant initial BH spin a,i n itiai = 0.01 and the K05 align- 
ment criterion (upper panels), (ii) also with K05 and con- 
sidering an initial BH spin proportional to the BH mass, 
ainitial = Mbh, initiai/10 7 Mq (middle panels) , and (iii) adopt- 
ing the same initial BH spin as in (i) but without considering 
K05 (lower panels). Vertical solid lines represent the average 
BH spin using the full sample for each morphological type. 
For reference, average values from model B are shown by 
the dashed lines. 

In all cases considered, elliptical galaxies host both 
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Figure 4. Normalized distributions of BH spin in four different 
BH mass bins obtained from model A using a constant initial BH 
spin a.j n j t j a i = 0.01 and the K05 alignment criterion. Solid vertical 
lines indicate average values of the BH spin calculated from the 
distributions shown in each panel. For reference, we also show the 
average values obtained from model A without K05 (dotted lines) 
and model B with K05 (dashed lines). The dashed histogram 
shows the distributions of model A with K05 and initially spinless 
BHs. 



rapidly and slowly-rotating BHs (the distributions show two 
peaks, at a~ and a~ 1), with a marked preference for 
high rotation; low values of spin are associated to low stel- 
lar mass galaxies (Msteliar i$ 10 10 Mq), with BH growth his- 
tories characterised by a typically low number of accretion 
episodes. Spiral galaxies, on the other hand, show a larger 
population of low spin BHs. Even though spirals also show 
a peak at a~ 1, the average values of BH spin in ellipticals, 
for both models A and B, are always higher than the aver- 
age BH spin in spirals. This trend does not depend on the 
morphology criterion, since threshold values rthrcsh > 0.7 
produce practically the same spin distributions (i.e. varia- 
tions in the average values are smaller than 20%). Moreover, 
even a threshold as small as rthresh = 0.5 produces the same 
trend of increasing spin values for BHs hosted by elliptical 
galaxies. Our results change only slightly when the model 
does not consider the K05 criterion (lower panels), produc- 
ing an almost negligible decrease of the average BH spin for 
both elliptical and spiral galaxies. However, the effect of not 
applying the criterion K05 has strong impact in model B, 
producing very low spin values in all populations (i.e. accre- 
tion episodes randomly add/subtract angular momentum to 
the BH). 

This morphology-related dichotomy also arises from the 
analysis of BH spin distributions as a function of the BH 
mass. Fig. [4] shows the resulting normalized distributions 
from model A with an initial BH spin fiinitiai = 0.01 and 
K05 alignments for four different BH mass ranges. The dis- 
tributions are characterised by a progressive increase of the 
relative number of BHs with high spin values as we move 
to higher BH masses. This is clearly shown by the average 



values of BH spins represented by the vertical solid lines. 
Taking into account that our mod el reproduces the well 
known BH-bulge mass relation (e.g.. | Magorrian et al ] |l998l . 
iGebhardt et alj|2000l . lHaring fc Rixll2004l ). the link between 
massive elliptical galaxies and high BH spins is also recov- 
ered (see Fig. |3}. The similar average values denoted by 
solid (model A with K05) and dotted (model A without 
K05) lines support our previous conclusion about the minor 
influence of the K05 alignment criterion on the BH spins. 
Although the average values from model B (random angles 
with K05; dashed lines) are different from those given by 
model A, the general behaviours of the average values are 
si milar in both models. This is opp osite to previous results 
bv lKing. Pringle fc Hofmannl l2008. where the most massive 
objects with BH spins obtained by considering random an- 
gles show the inverse trend. The distributions for model A 
with K05 and with a null initial BH spin are shown as dashed 
histograms and, as can be seen, there are little differences 
with respect to the non-zero spin case (solid histogram), and 
the increase in the average spin value as a function of BH 
mass is also obtained. 

The three different scenarios explored here imply that 
massive ellipticals ho s t mor e rapidly rotating BHs than spi- 
rals, as ISikora et al.l l|2007h proposed in their revised ver- 
sion of the "spin paradigm". This fact, together with the 
lack of a dependence ol the z = BH spin distributions 
on their initial value, confirm that the accretion and assem- 
bly history almost entirely define the final BH spin value. 
This was suggested earlier by our analysis of the BH growth 
tracks (Fig. [2}. This result is contrary to what is found 
from randomly oriented accretion discs, since in this case 
the final distribution retains the memory of its initial condi- 
tions, as has already been reported in previous works (e.g. 
IVolonteri 6^120051 ) . We also notice that the K05 criterion 
has a small impact on the final average value of the BH 
spin. In the following section we discuss into more detail the 
effects of considering K05 by using an individual example 
galaxy extracted from the model. 



4 BLACK HOLE SPIN DEVELOPMENT 

In this section we present in detail the assessment of the 
contribution from each BH growth mechanism to the final 
spin value. We also evaluate the influence of assumptions 
regarding the accretion disc orientation and the effect of 
K05 alignments. 

Fig. [5] shows the contribution of each BH growth mech- 
anism to the BH spin value at z = including mergers 
between BHs, and gas accretion driven by gas cooling pro- 
cesses, galaxy minor/major mergers, and disc instabilities in 
the host galaxy. The calculation of the contributions from 
each mechanism is explained in detail in Section[2]3] The left 
panel shows the results obtained from model B with a con- 
stant initial BH spin a.; n iti a i = 0.01 and the K05 alignment 
criterion. We can see that the main spin-up mechanisms are 
disc instabilities and galaxy major mergers, the same pro- 
cesses that provide the main source of BH growth (Fig. 2 
in LCP08). This result is not necessarily expected a priori 
since, in model B, the gas involved in minor and major merg- 
ers comes from different galaxies with random galaxy disc 
orientations, which may originate a misalignment of the ac- 
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Figure 5. Contribution to the final BH spin from the five different mechanisms involved in the BH growth: mergers between BHs 
(dot-dashed line), and gas accretion driven by galaxy major mergers (triple dot-dashed lines), galaxy minor mergers (dashed line), galaxy 
disc instabilities (solid line), and gas cooling (dotted lines). Left panel: results from model B with a constant initial BH spin a; n ; t ; a i = 0.01 
and the K05 alignment criterion. Middle panel: results from model A, with ai n i t i a i = 0.01 and K05. Right panel: same as middle panel 
but without the K05 criterion. 



cretion disc. Note that gas cooling processes produce a minor 
average contribution throughout the whole BH spin range, 
mainly due to the low mass accretion rates during this pro- 
cess in addition to the random relative orientations, which 
result in roughly alternating additions and subtractions to 
the BH spin modulus. 

The other two plots of Fig. show results from model 

A, with and without considering K05 alignments (middle 
and right panels, respectively). In both cases, we set a con- 
stant initial BH spin ainitiai — 0.01. We can see that, as 
in model B (left panel), disc instabilities are also one of 
the main spin-up mechanisms. However, the main difference 
with respect to the case of randomly oriented accretion discs 
is the important influence gained by gas cooling processes. 
The contribution of minor mergers also becomes more im- 
portant than in model B, being comparable to that of major 
mergers. 

These results can be understood recalling that, during 
galaxy mergers, the angular momentum of the accretion disc 
arises as a result of the contribution from the gas mass- 
weighted angular momenta of the merging galaxies, while 
in self-interaction processes, like gas cooling and disc in- 
stabilities, the accretion disc simply takes the direction of 
the host halo angular momentum or the galaxy disc, re- 
spectively. This indicates that, at high redshifts, the first 
non-aligned accretion discs are likely to come from galaxy 
mergers. Taking into account the comparable influence of 
disc instabilities and gas cooling processes on the results of 
model A, our conclusions on the development of the BH spin 
based on this model would not be strongly affected by the 
current modeling of disc instabilities, which is a simplified 
version of the actua l process taking place in real galaxies 
l| Athanassoulal [2008 ) . This is not necessarily true for model 

B, where disc instabilities are major contributors to the BH 
spin. 

Finally, mergers between BHs make a low relative con- 
tributions to the final value of the BH spin in all the 
three cases analysed. This is in agreement with previous 
works that indicate that high spins are achieved by accre- 
tion of ba ryonic matter rathe r than by mergers between 
BHs (e.g. [Hughes fc Blanfordl 120031 ; iPeirani fc de Freitasl 



20081; iBerti fc Volonteril |200S|; [Moderski fc Sikoral Il996l ; 



Moderski. Sikora fc Lasota 19981 ). We remark that the rel- 
ative importance of each growth mechanism is not affected 
by changes in the initial BH spin values tested in this work. 

These results are only slightly modified when K05 align- 
ments are not used, in which case the angular momentum of 
the disc is simply vectorially added to the BH spin without 
considering alignments (right panel of Fig. [5}. In this case, 
gas cooling processes and disc instabilities contribute to the 
BH spin in a rather similar way to the case where the K05 
criterion is considered. This indicates that the angular mo- 
mentum of accretion discs is mainly distributed within nar- 
row cones around the BH spin direction, implying a slightly 
increased final spin value for a model with K05. In the more 
rare case of counter-alignments, the accretion episodes can 
also produce more important spin-down effects. 

Regarding the effect of K05 alignments, we focus on an 
individual case in order to illustrate its effect on the BH spin 
development. Fig. [S] shows the evolution of several proper- 
ties of the BH-disc system of one example galaxy taken from 
model A, with and without considering the K05 alignment 
criterion (left and right panels, respectively). The upper pan- 
els show the relative contribution to the angular momentum 
of the BH, 5J/J BH (where S J = Jig? 1 - J^ tial is the an- 
gular momentum variation between consecutive snapshots 
of the underlying A^-Body simulation and Jbh = Jbh' 11 ^ 
coming from accretion discs formed by the different gas ac- 
cretion events (see Fig. [5}. In both cases, with and without 
K05, some accretion episodes contribute to the BH angu- 
lar momentum in amounts comparable or even higher than 
the current value of Jbh- The middle left panel shows the 
angle between the angular momentum of the accretion disc 
and the BH spin; the middle right panel presents the vari- 
ation in this angle arising from not including the K05 cri- 
terion. Lower panels depict the mass accreted by the BH 
in each snapshot (left-hand y-axis) and the BH spin his- 
tory (right-hand y-axis). Notice that the inclusion of K05 
induces important changes in the details of the BH spin de- 
velopment. In particular, the accretion episodes that take 
place between 2.5 ^ z ^ 3.5 do not change significantly 
the BH spin value when K05 is not applied, whereas when 
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Figure 6. Rcdshift evolution of properties of the BH-disc system for an example galaxy extracted from SAG (spiral galaxy with 2 = 
stellar mass 

^Stellar = 8 X 10 9 M Q and a BH mass of Mbh = 2 x 10 7 Mq ). The left and right panels show the results from model A 
with and without the K05 alignment criterion, respectively. The upper panels show the relative contribution of the angular momentum 
of the accretion discs to the angular momentum of the BH before the accretion, Jbh = Jgn' 3 ', giving a variation 5 J = Jg^ 1 — Jbh 
between consecutive snapshots of the iV-Body simulation. The middle left panel presents the angle between the accretion disc and the 
BH spin, 0; the middle right panel shows the variation in <j> when the K05 criterion is removed. Lower panels show the mass accreted by 
the BH and the BH spin development in solid and dashed lines, respectively. 



considering K05 the BH spin shows important variations. 
However, there is no final important net difference between 
these two models, since the sequence of accretion episodes 
that increase/decrease the BH spin almost cancel each other 
completely. Note that applying K05 also has a small effect 
on the resulting angles between BH spins and the galaxy an- 
gular momenta, producing differences of less than 4 degrees. 
Finally, it is interesting to notice, by looking at the middle 
panels, that in ~ 60 per cent of the accretion episodes the 
angle between the angular momenta of the disc and the BH 
is 4> < 90 degrees. 



Our inclusion of accretion disc orientations arising from 
the LSS, with or without K05 alignments, produces impor- 
tant changes in the contributions of different processes to the 
BH spin development, with respect to the results obtained 
using random orientations. This could have an important 
impact in our understanding of the connection between the 
BH spin and t he luminous activity of A GN. For instance, re- 
cent studies bv lCiotti fc Qstrikerl jjoO^ ) show that RL QSOs 
might be produced by recycled gas rather than mergers. In 
our model, these objects would be characterised by high spin 
values (see for instance Fig. [5]l ; however, if we considered a 
Monte-Carlo model, these objects would probably show a 
low spin (e.g. lKing fc Pringldl2007h . leading to a completely 
different interpretation of the physical phenomenon at work. 



5 BEHIND THE NATURE OF AGN RADIO 
LOUDNESS 

The results of our model have shown that the properties 
of black holes and their host galaxies are closely connected, 
in agreement with numerous observational and theoretical 
works on this subject. However, our model does not pro- 
vide direct information on the radio loudness of the galaxy 
population. In an attempt to classify a galaxy as a radio- 
loud object, we explore different possibilities using available 
properties of galaxies and their BHs provided by the model. 
Among the most important quantities required to fully char- 
acterise the BH and its accretion disc are the way in which 
gas accretion proceeds (extensively discussed in Section [1]), 
the BH mass, the BH spin, and the magnetic field. 

Observational results indicate that the average BH 
mass in radio-loud systems is about twice as la rge as 
radio-quiet 
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Metcalf fc Magliocch ctti 2006). Taking into account the well 
known relationship between BH and host bulge mass, it is 
natural to extend the dependence of radio loudness to the 
stellar mass of the host galaxy. Nonetheless, both popula- 
tions of RL and RQ QSOs have been detected in massive 
elliptical galaxies (e.g. lFlovd et al.ll2004h : also, BHs of equal 
mass show differing radio luminosities by up to ~ 4 orders o f 
magnitude l|Hol2002l . iDunlop et albooa . ISikora et all2007l ) . 
Therefore, it can be inferred that the influence of this BH 
property on the jet production, and the consequent strength 
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Figure 7. Fraction of z < 0.1 MRL galaxies as a function of stellar mass for model A (left panel) and B (right panel), both using an 
initial BH spin ^initial = 0.01 and including K05 alignments. Blue lines show the fractions for samples with Eddington ratios A < AmaXi 
with Amax = 10 — 1 , 10 — 2 , 10 — 3 and 10 — 4 (top to bottom lines); black lines, the fractions for samples with BH spin a > <2min> with 
Qmin = 0.3, 0.6 and 0.9 (top to bottom lines); red lines show the fractions resulting from combinations of limits on both accretion rate 
and BH spin. For reference, we show the fraction of active galaxies as a function of stellar mass (black dashed lines). 



of the radio emissio n, is secondary. Other results i n this di- 
rection are found bv lMerloni. Heinz fc di Matted (|2003l ). 

The second parameter that may play a role in the 
RL/RQ dichotomy is the amplitude of the local magnetic 
field from accretion flows which, in some cases, can produce 
large -scale fields that would help collimate strong relativistic 
jets ijKrolik fc Hawlevl 120071 ). It has been argued, however, 
that t he presence of such fields may not be of primary impor- 
tance |Spruit| [2008). In microquasars, that share the same 
physical principles than AGN on much smaller scales, it is 
possibl e that cold winds are responsible for the collimation 
of jets (|Bosch-Ramon. Romero fc Paredesll200d ). This alter- 
native further diminishes the potential relevance of magnetic 
fields as a primary collimation mechanism. 

The third physical property that can be linked to 
the radio loudness phenomenon is the BH spin as h as 
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2007). Fur- 



thermore, a recent study bv lLindner fc Fragile! (|2007l ) shows 
that the jet direction is influenced by the BH spin. There- 
fore, in addition to the accretion rate, we consider the BH 
spin as another possible candidate to break the degeneracy 
between the RQ and RL populations. 

With these caveats in mind, we assume that radio loud- 
ness is defined by the BH spin, the gas accretion rate ("ac- 
cretion paradigm" ) , or combin ations of both para meters (as 
in the "revised spin paradigm" . ISikora et al.l l 2007). We eval- 
uate the possibility that radio loudness is defined by these 
criteria using observational results on the fraction of radio- 
loud galaxies (/rl) w hich is though t to increase with both, 



BH and stellar mass (ICorbinlll997l; lLaorl Eoool; lLacv et all 
200ll;lBoro son 2002; Jarvis fc McLurell2002l ; lBest et al.ll2005l ; 
Metcalf fc Magliocchettil 120061 ). These relations are still a 
matter of debate, since it has been argued that they can 
simply be the resu lt of observational biases (see for instance 
I White et al.ll2007l ). In addition, since there are no measure- 
ments of /rl using complete samples of galaxies to date, we 



are limited to make qualitative comparisons to the observed 
data. The most complete st udy performed to date was pre- 
sented bv lBest et al.l (|2005h based on the Sloan Digital Sky 
Survey (SDSS) but, even in this case, they only consider FR 
Class 1 radio sources, characterised by low radio luminosi- 
ties. We will attempt to distinguish between the predictions 
of models A and B by using only qualitative comparisons. 

To this aim, we focus on z ~ active galaxies (i.e. with 
Eddington ratios A = LeH/^Edd > 0), and separate them 
into samples of model radio-loud galaxies (MLR) selected 
using different upper limits in accretion rates, A < Amax, 



with Amax = 10 



Kl- 



lO 3 and 10 4 , and three nor- 



malized minimum BH spin values, ayim = 0.3, 0.6 and 0.9. 
We test whether these MRL samples are able to mimic the 
observed radio-loud population. Fig. [7] shows the resulting 
fractions of MRL galaxies, /mrl, as a function of stellar mass 
for models A and B (cf. Section r2.4|l using an initial BH spin 
Sinitiai = 0.01 and including K05 alignments. Fractions ob- 
tained using thresholds in Eddington ratio only, A < Amax, 
are represented by blue lines. Black lines show the results for 
galaxies selected using a > Smin, and red lines, the results 
for samples obtained from the combination of limits on both 
spins and accretion rates. For reference, the results for the 
full population of active galaxies are represented by a black 
dashed line. 

As can be seen, for model A, the limits on Eddington 
ratios produce very similar fractions to what is obtained for 
the full AGN sample, mostly due to the lack of BHs with 
high A values at low redshift. The main features for samples 
limited by A are an increasing /mrl with stellar mass for 
M st eli ar < 10 9 A'/ Q , with a plateau up to Msteliar ~ 1O 1O M . 
Stringent limits, Amax < 10 -3 , produce a decreasing /mrl 
at the high stellar mass end. However these limits are very 
restrictive since observed RL galaxi es can show Eddington 
ratios as large as A w 10~ 2 (e.g. H3 120021 ; ISikora et al.l 
2007). This aspect should be taken into account in order 
to test these particular predictions. These trends are at 
odds with the monotonically increasing observational rela- 
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tion mentioned above. There is little effect on these trends 
when considering the case of a null initial BH spin. 

On the other hand, when considering thresholds on the 
BH spin alone, the fractions of MRL galaxies always increase 
with stellar mass. The results for the combined conditions, 
a > ayim and A < Amax, are only shown for Amax = 10 -1 
and 1CT 2 . These cases also show increasing relations between 
/mrl and stellar mass because of the strong influence of the 
limits on the BH spin. Finally, for all the selection criteria 
analysed, model B (right panel of Fig. [7}, characterised by 
random orientations, fails to produce monotonically increas- 
ing fractions of MLR galaxies as a function of stellar mass. 
When we consider null initial BH spins, the only appreciable 
change in the fractions of MLR galaxies is a slight increase 
in the slope of the relation with respect to the results shown 
in Fig. [3 

Fig. |8]shows the fractions of MRL galaxies as a function 
of BH mass for model A (lines colours are as in Fig. [7}. As 
can be seen, the fractions of galaxies selected according to 
the three sets of criteria on BH spin and accretion rates 
show similar behaviours to those obtained as a function of 
stellar mass. The same similarity between the dependences 
with BH and stellar mass arise for model B, with milder or 
negligible trends of the MRL fractions with BH mass. 

All the previous results correspond to the local Uni- 
verse. From our results on the development of BH spin (see 
for instance Fig.[SJ), it is natural to expect smaller fractions 
/mrl at higher redshifts, particularly when using the BH 
spin as the main parameter to define radio-loudness. How- 
ever, the fractions of active galaxies increase importantly at 
higher redshifts producing a combined effect. Galaxies host- 
ing BHs with masses Mbh > 10 s Mq lead to higher fractions 
/mrl with respect to those with BHs of the same mass in 
the local Universe, while galaxies with BHs masses in the 
range Mbh < 10 s M show the expected decreased frac- 
tions. As the redshift increases the relation becomes steeper 
as a function of both BH and stellar mass. This prediction 
can be confirmed by further observations of high redshift 
radio- sources. 

We have combined two models for the orientation of ac- 
cretion discs with three sets of criteria related to the accre- 
tion and spin paradigms, with the aim of isolating the main 
parameters that characterise the radio loudness of galaxies. 
Our results indicate that the dependence of the fractions 
of MRL galaxies with BH and stellar mass are a qualita- 
tive match to the increasing trend showed by observational 
results when thresholds on the BH spin are applied. It is 
important to notice that this is more naturally achieved in 
a model with LSS oriented accretion discs. The most impor- 
tant aspect to emphasize from our results is that this quali- 
tative agreement preferentially favours the "spin paradigm" 
as the likely scenario for radio loudness in the Universe. We 
remind the reader that the observed trend is still in debate. 
Therefore, further observational data on AGN are required 
in order to reach a proper understanding of the biases in- 
volved in the measurements of fractions of RL galaxies. Only 
once this problem is solved it will be possible to attempt to 
tip the scales in favour a single model, for example A or B. 
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Figure 8. Fraction of MRL galaxies as a function of BH mass 
for model A using an initial BH spin a; n ; t ; a i = 0.01 and including 
K05 alignments. Line types and galaxy samples are as in Fig. [7] 



6 SUMMARY AND PERSPECTIVES 

We have implemented a simple calculation of the BH spin 
in the semi-analytic model of galaxy formation described 
by LCP08 to study relationships between the BH spin, the 
BH mass and the host galaxy properties. The main advan- 
tage of using a semi-analytic model relies in the ability to 
follow, for the first time, the development of the BH spin 
from the contribution of different BH growth mechanisms, 
all immersed within a ACDM cosmology. Such processes are 
gas cooling, disc instabilities, galaxy mergers and mergers 
between BHs. Additionally, the model allows us to relate 
phenomena taking place at very small scales (BH and accre- 
tion disc) to those at larger scales (galaxy discs); we infer 
the orientation of the accretion disc using that of the source 
of accreted material (referred to as model A). Thus, the an- 
gular momentum of the host halo is used when the accretion 
is driven by gas cooling, while that of the galaxy disc is con- 
sidered during disc instability events; in the case of mergers, 
gas mass-weighted averages of the angular momentum of 
intervening galaxies are adopted. We also consider random 
accretion disc orientations (mod el B), as has been adopted 
in several previous works (e.g. IVolonteri. Sikora fc Lasotal 
12001 iKing fc Prin"gie]|2007f l . 

The main results and conclusions of this work are sum- 
marised as follows. 

• Massive bulges, and therefore giant elliptical galaxies, 
host almost exclusively rapidly rotating BHs. Although this 
effect is more pronounced in model A than in model B, 
in both models elliptical galaxies host, on average, BHs 
with higher spin values than spiral galaxies, i n accordance 
with t he revised "spin paradigm" proposed bv lSikora et al.l 
(2007). Therefore, the morphology-related bimodality of ra- 
dio loudness could be explained by the BH spin. 

• High mass BHs (Mbb. > 10 7 A/q) in model A, show 
little dependence on the initial BH spin assigned to them, 
indicating that their spin is almost entirely determined by 
the accretion history. Additionally, including the K05 align- 
ments of accretion discs only produces very mild changes 
in the final BH spin distribution. Results from model B are 
found to be more dependent on the initial values of the BH 
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spin and strongly dependent on whether the K05 alignment 
is applied or not. 

• For model A, gas cooling processes and disc instabilities 
are the main contributors to the final BH spin. It is impor- 
tant to bear in mind that gas cooling processes produce very 
low relative accretion onto the BH. On the other hand, in 
model B, the main mechanisms for BH spin-up are also those 
responsible for the highest accretion rates, namely, galaxy 
mergers and disc instabilities. This represents the main dif- 
ference between the two schemes adopted for the orientation 
of accretion discs. 

• We use measurements of the dependence of the fraction 
of RL galaxies as a function of stellar and BH masses in 
order to test different parameter combinations to select a 
population of model galaxies representative of the observed 
RL galaxy population. These parameters are assumed to 
be the accretion rate and/or the BH spin. We test three 
different criteria, using (i) an upper limit on the Eddington 
ratio ("accretion paradigm"), (ii) a lower limit on the BH 
spin, and (iii) simultaneous limits on both quantities (similar 
to the revised version of the "spin paradigm"). For model A, 
we find that criterion (i) predicts a roughly constant fraction 
of possible radio-loud galaxies as a function of BH and stellar 
mass (for M* > 1O 9 M0), while criterion (ii) produces a 
monotonically increasing relation. Finally, criterion (iii) also 
produces an increasing relation, although mostly due to the 
strong influence of the limit on the BH spin. For model B, we 
find that none of the proposed criteria are able to produce 
an increasing fraction of radio-loud galaxies with either the 
stellar or BH mass. 

• Taking into a ccount recent observational results (e.g . 
Best et al.l 120051; IChiaberge. Capetti fc Macchettol 120051; 



Metcalf fc Magliocchettf 2006 ) reporting an increasing frac- 
tion of radio-loud galaxies with BH and stellar masses, we 
find that the most suitable par adigm to explain r adio loud- 
ness is the "spin paradigm" l|Sikora et al.l [20071 ). Still, it 
should be taken into account that the se measurements may 
be subject to obse rvational biases fe.g. lCirasuolo et al 1l2ooa . 
IWhite et al.ll2007D . 

The implementation of a BH spin model presented here 
can be used as a new tool to understand the role of different 
processes on both very large and very small scales, within an 
ab-initio model of galaxy formation. Within our modeling we 
are able to reproduce results from Monte-Carlo attempts at 
describing BHs and the radio-loud galaxy population. Our 
analysis also singles out possible systematic biases in the 
model predictions resulting from the use of random orienta- 
tions for the accretion discs which continuously form around 
BHs. Our finding of a possible influence of the LSS align- 
ments on the detailed spin-up history of BHs can now also 
be considered in simple models, and be studied into more de- 
tail. However, new and improved observational constraints 
are still required to help distinguish between possible causes 
for the AGN radio activity in the Universe. 
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